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Abstract 

The production of jets is studied in deep-inelastic e^p scattering at large negative four 
momentum transfer squared 150 < < 15000 GeV^ using HERA data taken in 1999- 
2007, corresponding to an integrated luminosity of 395 pb~^. Inclusive jet, 2-jet and 3-jet 
cross sections, normalised to the neutral current deep-inelastic scattering cross sections, are 
measured as functions of Q^, jet transverse momentum and proton momentum fraction. 
The measurements are well described by perturbative QCD calculations at next-to-leading 
order corrected for hadronisation effects. The strong coupling as determined from these 
measurements is a^(Mz) = 0.1168 ± 0.0007 (exp.) tomto (th-) ± 0.0016 (PDF). 
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1 Introduction 



Jet production in neutral current (NC) deep-inelastic scattering (DIS) at HERA provides an 
important testing ground for Quantum Chromodynamics (QCD). While inclusive DIS gives 
only indirect information on the strong coupling via scaling violations of the proton structure 
functions, the production of jets allows a direct measurement of ag. The Bom level contribution 
to DIS (figure[I^) generates no transverse momentum in the Breit frame, where the virtual boson 
and the proton collide head on [1]. Significant transverse momentum Pt in the Breit frame is 
produced at leading order (LO) in the strong coupling by the QCD-Compton (figure [Ih) and 
boson-gluon fusion (figure [l]:) processes. 

In leading order the proton's momentum fraction carried by the emerging parton is given 
by ^ = XBj(l + M^2/Q^)- The variable xb^ denotes the Bjorken scaling variable, M12 the 
invariant mass of two jets of highest and the negative four momentum transfer squared. 
In the kinematical regions of low Q"^, low Pt and low ^, boson-gluon fusion dominates the jet 
production and provides direct sensitivity to the gluon component of proton density functions 
(PDFs) [2]. 




Figure 1: Deep-inelastic lepton-proton scattering at different orders in a^: (a) Bom contribution 
C(l), (b) example of the QCD Compton scattering C(a<j) and (c) boson-gluon fusion 0{as). 

Analyses of inclusive jet production in DIS at high Q"^ were previously performed by the 
HI [3] and ZEUS [4] collaborations at HERA. These analyses are based on data taken during 
1999 and 2000 (HERA-I) and use jet observables to test the mnning of the strong coupling and 
extract its value at the Z° boson mass. In this paper an integrated luminosity six times larger 
than available in the previous HI analysis [3] is used. The ratios of jet cross sections to the 
corresponding NC DIS cross sections, henceforth referred to as normalised jet cross sections, 
are measured. These ratios benefit from a partial cancellation of experimental and theoretical 
uncertainties. The measurements are compared with perturbative QCD (pQCD) predictions at 
next-to-leading order (NLO) corrected for hadronisation effects, and as is extracted from a fit of 
the predictions to the data. The measurements presented in this paper supersede the previously 
published normalised jet cross sections in [3]. 
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2 Experimental Method 



The data sample was collected with the HI detector at HERA in the years 1999 to 2007 when 
HERA collided electrons or positron^ of energy Ee = 27.6 GeV with protons of energy Ep = 
920 GeV, providing a centre-of-mass energy ^/s = 319 GeV. The data sample used in this 
analysis corresponds to an integrated luminosity of 395 pb^^, comprising 153 pb^^ recorded in 
e^p collisions and 242 pb^^ in e^p collisions. 

2.1 HI detector 

A detailed description of the HI detector can be found in [5, 6]. HI uses a right-handed co- 
ordinate system with the origin at the nominal interaction point and the z-axis along the beam 
direction. The positive z direction, also called the forward direction, is given by the outgoing 
proton beam. Polar angles 9 and azimuthal angles (p are defined with respect to this axis. The 
pseudorapidity is related to the polar angle 6* by 77 = — lntan(^^/2). The detector components 
important for this analysis are described below. 

The electromagnetic and hadronic energies are measured using the Liquid Argon (LAr) 
calorimeter in the polar angular range 4° < < 154° and with full azimuthal coverage 
[7]. The LAr calorimeter consists of an electromagnetic section (20 to 30 radiation lengths) 
with lead absorbers and a hadronic section with steel absorbers. The total depth of the LAr 
calorimeter varies between 4.5 and 8 hadronic interaction lengths. The energy resolution is 
as/E = 12%/ ^E / GeV® 1% for electrons and as/E = 50%/y/E / GeV ©2% for hadrons, 
as obtained from test beam measurements [8]. In the backward region (153° < 9 < 177°) en- 
ergy is measured by a lead/scintillating fibre Spaghetti-type Calorimeter (SpaCal) composed of 
an electromagnetic and a hadronic section [6]. The central tracking system (20° < 9 < 160°) 
is located inside the LAr calorimeter and consists of drift and proportional chambers, comple- 
mented by a silicon vertex detector covering the range 30° < 9 < 150° [9]. The chambers and 
calorimeters are surrounded by a superconducting solenoid providing a uniform field of 1.16 T 
inside the tracking volume. The luminosity is determined by measuring the event rate of the 
Bethe-Heitler process (ep — > ep'j), where the photon is detected in a calorimeter close to the 
beam pipe at z = —103 m. 

2.2 Event and jet selection 

The NC DIS events are triggered and selected by requiring a compact energy deposit in the 
electromagnetic part of the LAr calorimeter. The scattered electron is identified as the isolated 
cluster of highest transverse momentum [10]. Its reconstructed energy is requested to exceed 
1 1 GeV. Only the regions of the calorimeter where the trigger efficiency is greater than 98% 
are used for the detection of the scattered electron. These requirements ensure that the overall 
trigger efficiency reaches 99.5%. In the central region, 30° < 9 < 155°, the cluster has to be 
associated with a track measured in the inner tracking chambers and matched to the primary 

'Unless otherwise stated, the term "electron" is used in the following to refer to both electron and positron. 
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event vertex. The z-coordinate of the primary event vertex is required to be within ±35 cm of 
the nominal position of the interaction point. 

The remaining clusters in the calorimeters and charged tracks are attributed to the hadronic 
final state, which is reconstructed using an energy flow algorithm that avoids double counting of 
energy [11, 12]. Electromagnetic and hadronic energy calibration and the alignment of the HI 
detector are performed following the same procedure as in [10]. The total longitudinal energy 
balance, calculated as the difference of the total energy E and the longitudinal component of 
the total momentum P^, calculated from all detected particles including the scattered electron, 
must satisfy 35 < — < 65 GeV. This requirement reduces contributions of DIS events 
with hard initial state photon radiation. For the latter events, the undetected photons propa- 
gating in the negative z direction lead to values of this observable significantly lower than the 
expected value of 2Ee — 55.2 GeV. The E — requirement together with the scattered elec- 
tron selection also reduces contributions from photoproduction, estimated using Monte Carlo 
simulations. Cosmic muon and beam induced background is reduced to a negligible level after 
combining these cuts with the primary event vertex selection. Elastic QED Compton and lep- 
ton pair production processes are suppressed by rejecting events containing additional isolated 
electromagnetic deposits and low hadronic calorimeter activity. 

The kinematical range of this analysis is defined by 

150 < g2 < 15000 GeV^ and 0.2 < y < 0.7, 

where y = Q^/{sxBj) quantifies the inelasticity of the interaction. These two variables are 
reconstructed from the four momenta of the scattered electron and the hadronic final state par- 
ticles using the electron-sigma method [13]. The selection of events passing all the above cuts 
is the NC DIS sample, which forms the basis of the subsequent analysis. 

The jet finding is performed in the Breit frame, where the boost from the laboratory system 
is determined by Q^, y and by the azimuthal angle of the scattered electron. Particles of the 
hadronic final state are clustered into jets using the inclusive kr algorithm [14] with the massless 
Pt recombination scheme and with the distance parameter Rq = lin the i] — cf) plane. The cut 
—0.8 < r/Lab < 2.0, where r^Lab the jet pseudorapidity in the laboratory frame, ensures that 
jets are contained within the acceptance of the LAr calorimeter and are well calibrated. 

Jets are ordered by decreasing transverse momentum Pt in the Breit frame, which is iden- 
tical to the transverse energy Et for massless jets. The jet with highest Pt is referred to as 
the "leading jet". Every jet with the transverse momentum Pt in the Breit frame satisfying 
7 < Pt < 50 GeV contributes to the inclusive jet cross section. The upper cutoff is necessary 
for the integration of the NLO calculation. The steeply falling transverse momentum spectrum 
leaves almost no jets above 50 GeV. Events with at least two (three) jets with transverse mo- 
mentum 5 < Pt < 50 GeV are considered as 2-jet (3 -jet) events. In order to avoid regions of 
phase-space where fixed order perturbation theory is not reliable [15], 2-jet events are accepted 
only if the invariant mass Mu of the two leading jets exceeds 16 GeV. The same requirement, 
Mi2 > 16 GeV, is applied to the 3-jet events so that the 3-jet sample is a subset of the 2-jet 
sample. 

After this selection, the inclusive jet sample contains a total of 14381 1 jets in 104014 events. 
The 2-jet sample contains 47278 events and the 3-jet sample 7054 events. 
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2.3 Definition of the observables 



The measurements presented in this paper refer to the phase-space given in table [IJ Normalised 
inclusive jet cross sections are measured as functions of and double differentially as function 
of Q'^ and the transverse jet momentum in the Breit frame. Normalised 2-jet and 3-jet cross 
sections are presented as a function of . In addition the 2-jet cross sections are measured 
double differentially as function of and the average transverse momentum of the two leading 
jets (Pi-) = i ■ (P-^'^*^ + P^'^*^) or as function of Q"^ and of the proton momentum fraction ^. The 
3-jet cross section normalised to the 2-jet cross section as function of is also presented. 

The normalised jet cross sections are defined as the ratio of the differential inclusive jet, 2- 
jet and 3-jet cross sections to the differential NC DIS cross section in a given bin, multiplied 
by the respective bin width W in case of a double differential measurement as indicated by the 
following equations: 

ctnc daNc/d<3^ 

The normalised inclusive jet cross section can be viewed as the average jet multiplicity in a 
given Q'^ region and the normalised multi-jet cross sections as multi-jet event rates. 



2.4 Determination of normalised cross sections 

In each analysis bin the normalised jet cross section is determined as 

o, N 



C. (4) 



NC 



Here A^j denotes the number of inclusive jets or the number of 2-jet or 3-jet events, respec- 
tively, while A^NC represents the number of NC DIS events in that bin. The bin dependent 
correction factor C takes into account the limited detector acceptance and resolution. The cor- 
rection factors are determined from Monte Carlo simulations as the ratio of the normalised jet 
cross sections obtained from particles at the hadron level to the normalised jet cross sections 
calculated using reconstructed particles. 

The following LO Monte Carlo event generators are used for the correction procedure: 
DJANGOH [16], which uses the Color Dipole Model with QCD matrix element corrections as 
implemented in ARIADNE [17], and RAPGAP [18], based on QCD matrix elements matched 
with parton showers in leading log approximation. In both Monte Carlo generators the hadroni- 
sation is modelled with Lund string fragmentation [19]. All generated events are passed through 
a GEANT3 [20] based simulation of the HI apparatus and are reconstructed using the same pro- 
gram chain as for the data. Both RAPGAP and DJANGOH provide a good overall description 
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of the inclusive DIS sample. To further improve the agreement between Monte Carlo and data 
for the jet samples, the Monte Carlo events are weighted as a function of Q"^ and y and as func- 
tion of Pt and rj of the leading jet in the Breit frame. In addition, they are weighted as a function 
of Pt of the second and third jets when present [21]. After weighting, the simulations provide 
a good description of the shapes of all data distributions, some of which are shown in figure [2l 

The binnings in Q^, Pt and used to measure the jet observables are given in tabled The 
associated bin purities, defined as the fraction of the events reconstructed in a particular bin that 
originate from that bin on the generator level, are typically 70% and always greater than 60%. 
The correction factors deviate typically by less than 20% from unity, but reach 40% difference 
from unity in the bin 5 < {Pt) < 7 GeV for the 2-jet cross section. Arithmetic means of the 
correction factors determined from the reweighted RAPGAP and DJANGOH event samples are 
used and half of the difference is assigned as a model uncertainty. 

The above correction factors include QED radiation and electroweak effects. The effects 
of QED radiation, which are typically 5%, are corrected for by means of the HERACLES [22] 
program. The LEPTO event generator [23] is used to correct the e~^p and e^p data for their 
different electroweak effects which largely cancel in normalised jet cross sections leaving them 
below 3%. The resulting pure photon exchange cross sections obtained from e'^p and e~p data 
samples are then averaged. 

2.5 Experimental uncertainties 

The systematic uncertainties of the jet observables are determined by propagating the corre- 
sponding estimated measurement errors through the full analysis: 

• The relative uncertainty of the electron energy calibration is typically between 0.7% and 
1% for most of the events and increases up to 2% for electrons in the forward direction. 
The absolute uncertainty of the electron polar angle is 3 mrad. Uncertainties in the elec- 
tron reconstruction affect the event kinematics and thus the boost to the Breit frame. This 
in turn leads to a relative error of 0.5% to 1.5% on the normalised cross sections for each 
of the two sources, electron polar angle and energy. 

• The relative uncertainty on the energy of the total reconstructed hadronic final state as well 
as of jets is estimated to be 1.5% [21]. It is dominated by the uncertainty of the hadronic 
energy scale of the calorimeter. This error is estimated using a procedure similar to that 
used in [10] based on the transverse momentum conservation in the laboratory frame 
between the hadronic final state Pt^h and the electron e- This systematic uncertainty is 
reduced with respect to the previous measurement [3] due to the restricted pseudorapidity 
range in which jets are reconstructed and due to the improved statistics in the calibration 
procedure. The hadronic energy scale uncertainty affects mainly the jet cross section 
through the calibration of P^ and, to a lesser extent, the NC DIS cross section through 
the reconstruction of y. The resulting errors range between 1% and 5% and increase up to 
7% when P^ exceeds 30 GeV. The relative uncertainty due to the hadronic energy scale 
is reduced on average by about 20% for the normalised jet cross sections compared to the 
jet cross sections. 
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• The model dependence of the detector correction factors is estimated as described in 
section It reflects the sensitivity of the detector simulation to the details of the model, 
especially the parton showering, and their impact on the migration between adjacent bins 
in Pt. The model dependence ranges typically from 1% to 2% for Pt below 30 GeV and 
to 4% above, independently of Q"^. 

• The uncertainties of the luminosity measurements, the trigger efficiency and the electron 
identification efficiency cancel in the normalised cross section. In addition, the model de- 
pendence of the QED radiative corrections, which is estimated to be 1% [10], is expected 
to cancel in the normalised cross sections. 

The statistical errors for the normalised inclusive jet cross section take into account the 
statistical correlations which arise because there can be more than one jet per event [21]. The 
statistical errors are considerably smaller compared to the previous HERA-I publication [3]. 
They are typically between 1% and 2% for the normalised inclusive and 2-jet cross sections and 
do not exceed 10% in the regions of high transverse momentum Pt or high boson virtuality Q^. 

The dominant experimental errors on the jet cross sections arise from the uncertainty on 
the hadronic energy scale. The second most important source of systematic errors is the model 
dependence of the data correction, which becomes comparable to or exceeds the former in 
regions of highest jet Pt- The overall experimental error, calculated as the quadratic sum of 
all the contributions inventoried above, ranges typically between 3% and 6%, but increases up 
to 15% in the regions of highest Pt or Q^, dominated there by statistical uncertainties. The 
experimental errors for normalised cross sections are reduced by 30% up to 50% compared to 
those for unnormalised cross sections. 



3 NLO QCD prediction of jet cross sections 

Reliable quantitative predictions of jet cross sections in DIS require the perturbative calculations 
to be performed at least to next-to-leading order in the strong coupling. By using the inclusive 
kT jet algorithm with radius parameter R = 1, the observables used in the present analysis are 
infrared and coUinear safe and the non-perturbative effects are expected to be small [2]. In ad- 
dition, applying this algorithm in the Breit frame has the advantage that initial state singularities 
can be absorbed in the definition of the proton parton densities [24]. 

Jet cross sections are predicted at the parton level using the same jet definition as in the data 
analysis. The QCD predictions for the jet cross sections are calculated using the NLOJET-i-i- 
program at NLO in the strong coupling [25]. The NC DIS cross section is calculated at 0{as) 
with the DISENT package [26]. The FastNLO program [27] provides an efficient method to 
calculate these cross sections based on matrix elements from NLOJET-I-+ and DISENT, con- 
voluted with the PDFs of the proton and as a function of a^. The program includes a coherent 
treatment of the renormalisation and factorisation scale dependences of all ingredients to the 
cross section calculation, namely the matrix elements, the PDFs and a^. 

When comparing data and theory predictions the strong coupling at the Z'^ boson mass 
is taken to be as{Mz) = 0.1168 and is evolved as a function of the renormalisation scale 
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with two loop precision. The calculations are performed in the MS scheme for five massless 
quark flavours. The PDFs of the proton are taken from the CTEQ6.5M set [28]. The factori- 
sation scale is taken to be Q and the renormalisation scale to be a/ (<5^ + -P| obs)/2 
for the NLO predictions, with Pt,oU denoting Pt for the inclusive jet, (Pt) for 2-jet and 
I ■ (P^^*^ + P^^*^ + P^^^) for the 3-jet cross sections. This choice of the renormalisation scale 
is motivated by the presence of two hard scales, Pt and Q in the jet production in DIS. For the 
calculation of inclusive DIS cross sections, the renormalisation scale = Q is used. No QED 
radiation or exchange is included in the calculations, but the running of the electromagnetic 
coupling with Q"^ is taken into account. 

Hadronisation corrections are calculated for each bin using Monte Carlo event generators. 
These corrections are determined as the ratio of the cross section at the hadron level to the 
cross section at the parton level after parton showers. They typically differ by less than 10% 
from unity and are obtained using the event generators DJANGOH and RAPGAP which agree 
to within 2% to 4%. The arithmetic means of the two Monte Carlo hadronisation correction 
factors are used, while the full difference is considered as systematic error. 

DJANGOH and RAPGAP both use the Lund string model of hadronisation. The analytic 
calculations carried out in [29] provide an alternative method to estimate the effects of hadro- 
nisation and to cross-check the hadronisation correction procedure described above. They are 
based on soft gluon power corrections and result in a shift of the perturbatively calculated spec- 
trum of the inclusive jets: 

The size of the non-perturbative shift 5 (Pr)j^p can be calculated up to one single non-pertur- 
bative parameter «o(/^/) = /W/ ^ Jo' acs{k)dk, which is the first moment of the effective non- 
perturbative coupling Q;eff(A*) matched to the strong coupling as{fi) at the scale /x/. The value of 
ao{nj), expected to be universal [30], was measured to be ao{nj — 2 GeV) ^ 0.5 using event 
shapes observables in DIS by the HI Collaboration [31]. The hadronisation correction factors 
so calculated for the inclusive jet cross section differ in most of the bins by less than 2% from the 
average correction factor obtained from DJANGOH and RAPGAP and the maximum difference 
in all bins does not exceed 5% which is within the estimated uncertainty of the hadronisation 
correction. 

The dominant theoretical error is due to the uncertainty related to the neglected higher or- 
ders in the perturbative calculation. The accuracy of the NLO calculation is conventionally 
estimated by separately varying the chosen scales for // f and //j. by factors in the arbitrary range 
0.5 to 2. At high transverse momentum, above 30 GeV, the pQCD calculations do not depend 
monotonically on Hr in some bins. This happens in the two highest bins for the inclusive 
jet cross section and in six bins for the 2-jet cross section, where the largest deviation from 
the central value is found for factors well inside the range 0.5 to 2. In such cases the difference 
between maximum and minimum cross sections found in the variation interval is taken, in order 
not to underestimate the scale dependence. Renormalisation and factorisation scale uncertain- 
ties are added in quadrature, the former outweighing the latter by a factor of two on average. 
The uncertainties originating from the PDFs are estimated using the CTEQ6.5M set of parton 
densities. 
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Normalised jet cross sections are calculated by dividing the predicted jet cross sections by 
the NC DIS cross sections. The renormalisation scale uncertainties are assumed to be uncorre- 
cted between NC DIS and jet cross sections, as well as between 3-jet and 2-jet cross sections 
for their ratio, whereas the factorisation scale and the parameterisation uncertainty of the PDFs 
are assumed to be fully correlated. 



4 Results 

In the following, the normalised differential cross sections are presented for inclusive jet, 2-jet 
and 3-jet production at the hadron level. Tables [3] to [6] and figures [3] to [6] present the measured 
observables together with their experimental uncertainties and hadronisation correction factors 
applied to the NLO predictions. These measurements are subsequently used to extract the strong 
coupling as as shown in the table [9] and figures 171 to [T2l 

4.1 Cross section measurements compared to NLO predictions 

The normalised inclusive jet cross sections as a function of are shown in figure[3^ and table [3] 
together with the NLO predictions and previous measurements by HI based on HERA-I data 
[3]. For comparison, the HERA-I data points were corrected for the phase space difference due 
to the slightly smaller jet pseudorapidity range of the present analysis. The double differential 
results as a function of Pt in six ranges of are given in figure|4]and tablelH Normalised 2-jet 
(3-jet) cross sections as a function of Q"^ and their comparison to NLO are also shown in figure 
[3)3 (I3t) and table[3l while the ratio 3-jet to 2-jet is shown in figure[3]l. Figures [51 [6] and tables [51 
[6lpresent the normalised 2-jet cross section as a function of {Pt) and ^ in six ranges of Q^. 

The new measurement of the normalised inclusive jet cross section is compatible with the 
previous HI data. The precision is improved by typically a factor of two, as can be seen for 
example in figure [3^. The QCD NLO predictions for all normalised jet cross sections provide 
a good description of the data over the whole phase space. In almost all bins the theory er- 
ror, dominated by the /i, scale uncertainty, is significantly larger than the total experimental 
uncertainty, which is dominated by the hadronic energy scale uncertainty. 

The normalised inclusive jet cross section, which may be interpreted as the average jet 
multiplicity produced in NC DIS, increases with as the available phase space opens (figure 
[3^) as do the 2-jet and 3-jet rates (figure[3j5 and[3t). As Q'^ increases, the Pt jet spectra become 
harder as can be seen in figure [H and \5\ The 3-jet rate is observed to be nearly seven times 
smaller than the 2-jet rate as shown in figure [3li. The 2-jet rates measured as a function of 
and the momentum fraction ^ are well described by the NLO calculations (figure [6l). Kinematic 
constraints from the considered y range and the restricted invariant mass of the jets lead to a 
reduction of the 2-jet rate at low ^ and a rise at large ^ with increasing Q"^. 
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4.2 Extraction of the strong coupling 



The QCD predictions for jet production depend on as and on the parton density functions of the 
proton. The strong coupling as is determined from the measured normalised jet cross sections 
using the parton density functions from global analyses, which include inclusive deep-inelastic 
scattering and other data. The determination is performed from individual observables and also 
from their combination. 

QCD predictions of the jet cross sections are calculated as a function of ^^(/ir) with the 
FastNLO package using the CTEQ6.5M proton PDFs and applying the hadronisation correc- 
tions as described in section [3l Measurements and theory predictions are used to calculate a 
X^{ois) with the Hessian method [33], where parameters representing systematic shifts of de- 
tector related observables are left free in the fit. The shifts in the electron energy scale, electron 
polar angle and the hadronic final state energy scale found by the fit are consistent with the 
a priori estimated uncertainties. This method takes into account correlations of experimental 
uncertainties and has also been used in global data analyses [33, 34] and in previous HI publi- 
cations [3, 35]. The experimental uncertainty of as is defined by the change in as which gives 
an increase in of one unit with respect to the minimal value. 

The correlations of the experimental uncertainties between data points were estimated using 
Monte Carlo simulations: 

• The statistical correlations between different observables using the same events are taken 
into account via the correlation matrix given in tables |7] and [8l 

• It is estimated that the uncertainty of the LAr hadronic energy scale is equally shared 
between correlated and uncorrelated contributions [3,21], while that from the electron 
energy scale is estimated to be 3/4 uncorrelated [10]. 

• The measurement of the electron polar angle is assumed to be fully correlated [10]. 

• The model dependence of the experimental correction factors is considered as fully un- 
correlated after the averaging procedure described in section [231 

The sharing of correlated and uncorrelated contributions between the different sources of un- 
certainty has the following impact on the as determination: when going from uncorrelated to 
fully correlated error for each source, the fitted value of as typically varies by half the total 
experimental error and the estimated uncertainty by less than 0.1% of a^. 

The theory error is estimated by the so called ojfset method as the difference between the 
value of as from the nominal fit to the value when the fit is repeated with independent variations 
of different sources of theoretical uncertainties as described in section [3l The resulting uncer- 
tainties due to the different sources are summed in quadrature. The up (or down) variations are 
applied simultaneously to all bins in the fit. The impact of hadronisation corrections on is 
between 0.4% and 1.0%, while that of the factorisation scale amounts to 0.5%. The sensitivity 
of as to the renormalisation scale variation of the inclusive NC DIS cross section alone is typ- 
ically 0.5%. The largest uncertainty, of typically 3% to 4%, corresponds to the accuracy of the 
NLO approximation to the jet cross sections estimated by varying the renormalisation scale as 
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described in section [3l An alternative method to estimate the impact of missing orders, called 
the band method, developed by Jones et al. [36] was tested and, for the present measurement, it 
leads to a smaller uncertainty on of typically 2%. 

The uncertainty due to PDFs is estimated by propagating the CTEQ6.5M errors. The typical 
size of the resulting error is 1.5% for determined from the normalised inclusive jet or 2-jet 
cross sections and 0.8% when measured with the normalised 3-jet cross sections. This uncer- 
tainty is twice as large as that estimated with the uncertainties given for the MSTW2008nlo90cl 
set [37] which in turn exceeds the difference between values extracted with the central sets of 
CTEQ6.5M and MSTW2008nlo. The PDFs also depend on the value of a^. Potential biases on 
the as extraction from that source have been studied in detail previously [3]. For this analysis, 
the resulting uncertainty is found to be negligible. 

Individual fits of as{fir) are made to each of the 24 measurements of the normalised double 
differential inclusive jet cross section, as shown in figure |7^. These individual determinations 
show the expected scale dependence. Equivalently, the values at each scale /i,. can be related 
to the value of the strong coupling as{Mz) at the Z° mass as shown in figured 

Then as{Mz) is determined by a common fit to the normalised inclusive jet cross section 
in four Pt bins for each region in Q^. The resulting six values are evolved from the scale Mz 
to the average Q in that region (figure [TOk). Finally, a central value as{Mz) is extracted from 
a common fit to all 24 measurements and given in table |9l The result of evolving this value 
together with its associated uncertainty is also shown as the curve and surrounding band in 
figure [TOk. 

The same fit procedure of successive combination steps is applied to the 24 points of the 
normalised 2-jet cross section with (Pt) > 7 GeV (figure Uh, M and [TOb). The bins with 
5 < (Pt) < 7 GeV are not used for the extraction of the strong coupling since the theory 
uncertainty is significantly larger than in the other bins (figure [5]). The fit procedure is also 
applied to the 6 points of the normalised 3-jet cross section (figure [TOb). The normalised 3-jet 
cross section (figure[3l:), which is 0{a'^), is preferred to the ratio of the 3-jet cross section to the 
2-jet cross section (figure [Hi), which is 0{al), due to better sensitivity to the strong coupling. 
The three values of as{Mz) determined from the normalised inclusive jet (24 points), 2-jet (24 
points) and 3-jet (6 points) cross sections are given in table |9] with experimental and theoretical 
uncertainties. All obtained values are compatible with each other within two standard deviations 
of the experimental uncertainty. 

The impact of the choice of renormalisation scale on the central value of as{Mz) is studied 
in the case of the normalised inclusive jet cross section by repeating the fit procedure with 
/ir = Pt and jj,r = Q instead of fir = a/ {Q"^ + Pj)/'^- I^i the first case the central value of the 
as{Mz) is found to be approximatively 0.7% smaller and in the latter approximatively 1.5% 
bigger with respect to the nominal fit, a difference which is well inside the estimated theoretical 
uncertainties. Similar deviations are observed for the normalised 2-jet and 3-jet cross sections 
when = <5 is used instead of /x^ = \/ (Q^ + Pt obs) /2- To get information on the description 
of the data by the NLO calculations as a function of the renormalisation scale, the of the fit 
is studied in the case of the normalised inclusive jet cross section for different values of the 
parameter x,., defined by jj,r = Xr ■ \/ (Q^ + -Rr)/2. The results are shown in figure [HI where 
the OLs fit is repeated for different choices of Xr and the corresponding -y^ values are shown. 
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The lowest value is obtained for — 0.5 while Xr choices above 2.0 and below 0.3 are 
disfavoured. 

The sensitivity of the strong coupling determination procedure to the choice of the jet def- 
inition is tested for the normalised inclusive jet and 2-jet cross sections by repeating all the 
extraction procedure using the anti-fc^ metric [38] instead of kr, but keeping the recombination 
scheme and the distance parameter unchanged. The resulting central value of as{Mz) differs in 
both cases by less than 0.6% from the central value extracted using the kr metric. 

In each Q"^ region the values of as{Mz) from different observables are combined taking into 
account statistical and systematic correlations. The resulting values, evolved from the scale Mz 
to the average /i,. of the measurements in each Q'^ region, are shown in figure[T2l This visualises 
the running of ag for scales between 10 and 100 GeV and the corresponding experimental and 
theory uncertainties. All 54 data points are used in a common fit of the strong coupling taking 
the correlations into account with a fit quality x^/ndf = 65.0/53 (see table HI), which is also 
shown in figure [T2l 

The values of as{Mz) obtained in this way are also consistent with the world averages 
a,(Mz) = 0.1176 ± 0.0020 [39] and a,(M^) = 0.1189 ± 0.0010 [40], and with the previous 
HI and ZEUS determinations from inclusive jet production measurements [3,4] and multijet 
production [41]. The experimental error on as{Mz) measured with each observable typically 
amounts to 1%. The combination of different observables, even though partially correlated, 
gives rise to additional constraints on the strong coupling and leads to an improved experi- 
mental uncertainty of 0.6%. The experimental error on is independent of the choice of 
renormalisation scale within the variation used to determine the theoretical uncertainty. The 
total error is strongly dominated by the theoretical uncertainty due to missing higher orders in 
the perturbative calculation which is about 4%. 

5 Conclusion 

Measurements of the normalised inclusive, 2-jet and 3-jet cross sections in the Breit frame in 
deep-inelastic electron-proton scattering in the range 150 < < 15000 GeV^ and 0.2 < y < 
0.7 using the HI data taken in years 1999 to 2007 are presented. Calculations at NLO QCD, 
corrected for hadronisation effects, provide a good description of the single and double differ- 
ential cross sections as functions of the jet transverse momentum Pt, the boson virtuality as 
well as of the proton momentum fraction ^. The strong coupling as is determined from a fit of 
the NLO prediction to the measured normalised jet cross sections. The normalisation leads to 
cancellations of systematic effects, resulting in improved experimental and PDF uncertainties. 
The experimentally most precise determination of as{Mz) is derived from a common fit to the 
normalised jet cross sections: 

as{Mz) = 0.1168 ± 0.0007 (exp.) (th.) ± 0.0016 (pdf) . 

The dominating source of the uncertainty is due to the renormalisation scale dependence, which 
is used to estimate the effect of missing higher orders beyond NLO in the pQCD prediction. This 
measurement improves the experimental precision on determinations from other recent jet 
measurements at HERA [3,4]. The result is competitive with those from e+e^ data [40,42] and 
is in good agreement with the world average [39,40]. 
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NC DIS Selection 


150 < g2 < 15000 GeV^ 0.2 < y < 0.7 


Inclusive jet 


7 <Pt< 50GeV 




-0.8 < r^Eb < 2.0 


2-jet 


5 < pjf", pI''^ < 50GeV 


Mi2 > 16GeV 


3-jet 


5 < pjf", pjf*^ pjf*' < 50GeV 



Table 1: Selection criteria for the NC DIS and jet samples. 



bin number corresponding Q'^ range 



1 


150 < Q2 


< 


200 GeV^ 


2 


200 < Q2 


< 


270 GeV^ 


3 


270 < Q2 


< 


400 GeV^ 


4 


400 < Q2 


< 


700 GeV^ 


5 


700 < Q2 


< 


5000 GeV^ 


6 


5000 < Q2 


< 15000 GeV^ 



bin letter corresponding Pt or {Pt) range bin letter corresponding ^ range 



a' 5 < Pt < 7 GeV 

a 7 < Pt < 11 GeV A 0.006 < ^ < 0.020 

b 11 < Pt < 18 GeV B 0.020 < ,^ < 0.040 

c 18 < Pt < 30 GeV C 0.040 < ^ < 0.080 

d 30 < Pt < 50 GeV D 0.080 < ,^ < 0.316 



Table 2: Nomenclature for the bins in Q^, Pt for the inclusive jet or (Pr) for 2-jets and { used 
in the following tables. In case of the normalised 2-jet cross section, the bin a' in (Pt) is not 
used for the a. extraction. 
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Normalised inclusive jet cross section in bins of 











total 


total 


single contributions to correlated uncertainty 


tiadronisation 


hadronisation 


bin 


normalised 


statistical 


total 


uncorrelated 


correlated 


electron 


electron 


hadronic 


correction 


correction 




cross 


uncert. 


uncert. 


uncertainty 


uncert. 


energy scale 


polar angle 


energy scale 


factor 


uncertainty 




section 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 




{%) 


1 


2.39 10^^ 


0.7 


3.2 


2.4 


2.1 


0.6 


0.8 


1.9 


0.95 


0.6 


2 


2.69 10^1 


0.7 


3.0 


2.2 


1.9 


0.6 


0.6 


1.7 


0.94 


0.6 


3 


3.11 10^1 


0.8 


2.9 


2.3 


1.8 


0.6 


0.4 


1.6 


0.94 


0.8 


4 


3.62 10^1 


0.8 


2.7 


2.2 


1.7 


0.7 


0.3 


1.5 


0.94 


0.6 


5 


4.26 10^1 


0.9 


2.7 


2.3 


1.6 


0.9 


0.1 


1.3 


0.93 


1.7 


6 


5.02 10"^ 


3.2 


5.7 


5.2 


2.4 


2.2 


0.3 


0.8 


0.93 


3.0 


Normalised 2-jet cross section in bins of 


1 


8.81 10^^ 


1.0 


2.9 


2.2 


1.9 


0.4 


0.7 


1.7 


0.94 


1.1 


2 


1.01 10~^ 


1.1 


2.6 


2.1 


1.7 


0.3 


0.6 


1.5 


0.93 


1.3 


3 


1.19 10^^ 


1.1 


2.4 


1.9 


1.5 


0.2 


0.4 


1.4 


0.93 


1.3 


4 


1.41 10^1 


1.2 


2.5 


2.0 


1.4 


0.3 


0.4 


1.3 


0.91 


1.1 


5 


1.75 10^1 


1.2 


2.4 


2.1 


1.2 


0.1 


0.2 


1.1 


0.91 


2.9 


6 


1.97 10^1 


4.4 


7.7 


7.6 


1.1 


0.3 


0.2 


1.0 


0.91 


2.9 


Normalised 3-jet cross section in bins of 


1 


1.19 10^^ 


2.6 


5.1 


4.1 


3.1 


0.4 


1.3 


2.8 


0.85 


2.4 


2 


1.29 10^^ 


2.8 


5.1 


4.2 


2.9 


0.3 


0.7 


2.8 


0.84 


1.7 


3 


1.68 10^^ 


2.7 


4.6 


3.8 


2.6 


0.1 


0.9 


2.5 


0.83 


1.0 


4 


2.06 10^^ 


2.9 


4.7 


4.0 


2.5 


0.3 


0.8 


2.4 


0.82 


0.6 


5 


2.36 10^^ 


2.8 


6.6 


6.2 


2.3 


0.4 


0.4 


2.2 


0.81 


1.2 


6 


2.82 10"^ 


9.2 


18.7 


18.5 


2.3 


0.4 


0.3 


2.3 


0.75 


3.6 


3-jet cross section normalised to 2-jet cross section in bins of 


1 


1.36 10"^ 


2.7 


4.4 


4.3 


1.2 


0.2 


0.5 


1.1 


0.91 


1.5 


2 


1.28 10^1 


3.0 


4.7 


4.5 


1.4 


0.5 


0.1 


1.3 


0.90 


1.0 


3 


1.41 10^1 


2.9 


4.5 


4.3 


1.1 


0.2 


0.1 


1.1 


0.90 


0.9 


4 


1.46 10^1 


3.1 


4.8 


4.6 


1.2 


0.7 


0.3 


1.0 


0.90 


0.6 


5 


1.35 10^^ 


3.0 


5.1 


5.0 


1.2 


0.4 


0.2 


1.1 


0.89 


1.5 


6 


1.43 10^1 


9.8 


14.2 


14.1 


1.3 


0.3 


0.2 


1.2 


0.82 


3.3 



Table 3: Normalised inclusive jet, 2-jet and 3-jet cross sections in NC DIS measured as a func- 
tion of Q'^ . The measurements refer to the phase-space defined in table [IJ In columns 3 to 9 
are shown the statistical uncertainty, the total experimental uncertainty, the total uncorrelated 
uncertainty including the statistical one and the total correlated uncertainty calculated as the 
quadratic sum of the following three components: the electron energy scale, the electron po- 
lar angle uncertainty and the hadron energy scale uncertainty. The sharing of the uncertainties 
between correlated and uncorrelated sources is described in detail in section |4.2[ The hadro- 
nisation correction factors applied to the NLO predictions and their uncertainties are shown in 
columns 10 and 11. The bin nomenclature of column 1 is defined in tabled 
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Normalised inclusive jet cross section in bins of and Pt 











total 


total 


single contributions to correlated uncertainty 


tiadronisation 


tiadronisation 


bin 


normalised 


statistical 


total 


uncorrelated 


correlated 


electron 


electron 


hadronic 


correction 


correction 




cross 


uncert. 


uncert. 


uncertainty 


uncert. 


energy scale 


polar angle 


energy scale 


factor 


uncertainty 




section 


(%) 


(%) 


(%) 


(%) 


(%) 


{%) 


(%) 




(%) 


1 a 


1.53 10^^ 


0.8 


2.7 


2.1 


1.7 


0.6 


0.6 


1.4 


0.94 


0.7 


1 b 


6.93 10-^ 


1.2 


4.5 


3.5 


2.9 


0.6 


1.1 


2.6 


0.97 


0.3 


1 c 


1.53 10"^ 


2.5 


6.1 


4.7 


3.9 


0.6 


1.6 


3.5 


0.96 


0.6 


1 d 


1.93 10"^ 


7.2 


10.6 


9.7 


4.4 


0.2 


1.3 


4.2 


0.95 


1.8 


2a 


1.66 10~^ 


0.9 


2.6 


2.2 


1.4 


0.7 


0.4 


1.2 


0.93 


0.6 


2b 


8.10 10"^ 


1.3 


3.9 


3.0 


2.5 


0.6 


0.8 


2.3 


0.97 


0.4 


2 c 


1.97 10"^ 


2.6 


5.8 


4.6 


3.6 


0.3 


0.9 


3.5 


0.96 


0.9 


2d 


2.67 10"^ 


7.1 


10.2 


9.1 


4.6 


0.4 


0.5 


4.5 


0.97 


3.2 


3a 


1.82 10^^ 


1.0 


2.8 


2.4 


1.4 


0.7 


0.4 


1.1 


0.92 


0.7 


3b 


9.82 10"^ 


1.3 


3.5 


2.7 


2.2 


0.5 


0.4 


2.1 


0.97 


1.0 


3 c 


2.76 10"^ 


2.4 


5.5 


4.4 


3.3 


0.3 


0.8 


3.2 


0.96 


0.4 


3d 


3.11 10^-' 


7.0 


9.8 


8.5 


4.8 


0.1 


1.9 


4.4 


0.95 


3.2 


4a 


2.02 10"^ 


1.1 


2.4 


2.2 


1.1 


0.7 


0.2 


0.9 


0.92 


0.5 


4b 


1.16 10^^ 


1.3 


3.4 


2.8 


2.0 


0.8 


0.4 


1.8 


0.96 


0.5 


4c 


3.83 10"^ 


2.3 


5.9 


4.9 


3.3 


0.5 


0.7 


3.1 


0.97 


1.5 


4d 


5.28 10"-' 


6.3 


8.9 


7.9 


4.1 


0.3 


0.6 


4.1 


0.96 


2.7 


5a 


2.13 10"^ 


1.2 


2.4 


2.1 


1.1 


0.8 


0.1 


0.7 


0.90 


2.4 


5b 


1.42 10^1 


1.3 


3.3 


2.8 


1.7 


0.9 


0.1 


1.5 


0.96 


1.1 


5c 


5.91 10^^ 


2.0 


4.7 


3.8 


2.7 


0.9 


0.1 


2.6 


0.97 


0.3 


5d 


1.09 10~^ 


4.4 


7.4 


6.3 


3.9 


0.1 


0.3 


3.8 


0.96 


3.1 


6a 


2.32 10~^ 


4.3 


8.1 


7.8 


2.3 


2.2 


0.4 


0.4 


0.90 


3.9 


6b 


1.62 10"^ 


4.8 


7.5 


6.9 


2.9 


2.8 


0.4 


0.7 


0.94 


2.5 


6c 


8.14 10"^ 


6.7 


9.8 


9.4 


2.6 


2.2 


0.1 


1.4 


0.96 


0.8 


6d 


2.66 10"^ 


9.7 


19.0 


18.8 


3.1 


0.6 


0.5 


3.0 


0.97 


3.6 



Table 4: Normalised inclusive jet cross sections as a function of and together with their 
relative errors and hadronisation correction factors. Other details are given in the caption to 
table [3l The bin nomenclature is defined in tabled 
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Normalised 2-jet cross section in bins of and (Pt) 











total 


total 


single contributions to correlated uncertainty 


tiadronisation 


tiadronisation 


bin 


normalised 


statistical 


total 


uncorrelated 


correlated 


electron 


electron 


tiadronic 


correction 


correction 




cross 


uncert. 


uncert. 


uncertainty 


uncert. 


energy scale 


polar angle 


energy scale 


factor 


uncertainty 




section 


{%) 


(%) 


(%) 


{%) 


(%) 


(%) 


{%) 




(%) 


1 a' 


9.14 10~3 


3.2 


3.5 


3.5 


0.6 


0.5 


0.1 


0.3 


0.83 


2.5 


1 a 


4.40 10~^ 


1.4 


2.9 


2.6 


1.4 


0.5 


0.5 


1.2 


0.94 


1.4 


1 b 


2.77 10^^ 


1.7 


4.4 


3.6 


2.6 


0.4 


1.0 


2.3 


0.96 


1.4 


1 c 


6.28 10^-'* 


3.5 


6.8 


5.3 


4.2 


0.3 


1.9 


3.8 


0.96 


1.4 


1 d 


6.87 lO"* 


10.5 


12.6 


11.7 


4.5 


0.1 


0.8 


4.5 


0.95 


1.8 


2 a' 


1.04 10~^ 


3.5 


4.2 


4.1 


0.9 


0.2 


0.8 


0.2 


0.83 


1.1 


2a 


4.91 10~^ 


1.5 


2.8 


2.5 


1.2 


0.4 


0.5 


1.0 


0.94 


1.3 


2b 


3.26 10^^ 


1.8 


3.8 


3.0 


2.2 


0.4 


0.8 


2.1 


0.96 


1.9 


2c 


7.80 10"^ 


3.6 


6.3 


5.1 


3.6 


0.2 


1.1 


3.4 


0.96 


1.7 


2d 


1.05 10^^ 


10.1 


12.5 


11.5 


4.9 


0.4 


0.5 


4.8 


0.92 


3.3 


3 a' 


1.13 10"^ 


3.7 


3.9 


3.9 


0.3 


0.3 


0.1 


0.2 


0.80 


1.4 


3a 


5.56 10"^ 


1.6 


2.6 


2.4 


0.9 


0.3 


0.3 


0.8 


0.92 


0.3 


3b 


3.99 10"^ 


1.8 


3.4 


2.9 


1.9 


0.4 


0.3 


1.8 


0.97 


2.2 


3c 


1.10 10^^ 


3.3 


5.9 


4.8 


3.3 


0.1 


0.8 


3.2 


0.96 


1.0 


3d 


1.16 10"-'' 


10.1 


13.2 


11.8 


5.9 


0.2 


2.9 


5.2 


0.94 


2.7 


4 a' 


1.41 10^^ 


3.9 


4.1 


4.1 


0.2 


0.1 


0.1 


0.2 


0.79 


2.7 


4a 


6.13 10"^ 


1.8 


2.7 


2.6 


0.8 


0.3 


0.1 


0.7 


0.90 


0.1 


4b 


4.80 10^^ 


1.9 


3.5 


2.9 


1.8 


0.6 


0.5 


1.7 


0.96 


1.3 


4c 


1.57 10^^ 


3.2 


6.3 


5.6 


3.0 


0.2 


0.6 


2.9 


0.97 


1.2 


4d 


2.09 lO"-'' 


9.1 


12.7 


11.8 


4.7 


0.1 


0.7 


4.7 


0.96 


2.6 


5 a' 


1.53 10^^ 


4.2 


9.2 


9.2 


0.8 


0.7 


0.2 


0.0 


0.77 


2.3 


5a 


6.95 10"^ 


1.9 


2.6 


2.5 


0.5 


0.1 


0.2 


0.5 


0.89 


3.6 


5b 


5.98 10"^ 


1.9 


2.9 


2.6 


1.3 


0.3 


0.1 


1.3 


0.94 


1.1 


5c 


2.49 10~^ 


2.8 


4.7 


4.0 


2.5 


0.1 


0.5 


2.4 


0.97 


2.2 


5d 


4.34 10~^ 


6.4 


9.1 


8.1 


4.2 


0.4 


0.7 


4.1 


0.96 


2.5 


6 a' 


1.30 10~^ 


16.2 


36.2 


36.1 


2.4 


2.3 


0.6 


0.6 


0.73 


16.2 


6a 


7.47 10"^ 


7.1 


10.6 


10.6 


0.5 


0.4 


0.1 


0.3 


0.88 


1.0 


6b 


6.42 10^^ 


6.9 


8.0 


7.8 


1.5 


0.9 


0.3 


1.1 


0.93 


4.4 


6c 


3.36 10"^ 


9.2 


9.7 


9.5 


1.6 


0.7 


0.4 


1.3 


0.95 


1.8 


6d 


1.03 10^^ 


15.7 


19.2 


18.8 


3.6 


0.6 


0.5 


3.5 


0.97 


4.3 



Table 5: Normalised 2-jet cross sections as a function of and (Pt) together with their relative 
errors and hadronisation correction factors. Other details are given in the caption to table[3l The 
bin nomenclature is defined in table [2l 
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Normalised 2-jet cross section in bins of and i 











total 


total 


single contributions to correlated uncertainty 


tiadronisation 


hadronisation 


bin 


normalised 


statistical 


total 


uncorrelated 


correlated 


electron 


electron 


hiadronic 


correction 


correction 




cross 


uncert. 


uncert. 


uncertainty 


uncert. 


energy scale 


polar angle 


energy scale 


factor 


uncertainty 




section 


{%) 


(%) 


(%) 


{%) 


(%) 


{%) 


(%) 




(%) 


1 A 


4.3C 10"^ 


1.4 


2.5 


2.1 


1.3 


0.7 


0.5 


1.0 


0.95 


0.7 


1 B 


3.37 10"^ 


1.6 


3.6 


2.9 


2.1 


0.4 


0.8 


1.9 


0.93 


1.7 


1 C 


9.22 10"-' 


2.9 


6.1 


4.8 


3.8 


0.4 


1.6 


3.4 


0.92 


3.2 


1 D 


1.88 10"-' 


6.5 


10.5 


8.8 


5.6 


1.2 


2.4 


4.9 


0.91 


1.3 


2 A 


4.20 10"^ 


1.7 


2.8 


2.5 


1.3 


0.9 


0.6 


0.7 


0.95 


0.6 


2 B 


4.44 10"^ 


1.6 


2.8 


2.3 


1.7 


0.2 


0.6 


1.6 


0.93 


2.7 


2C 


1.26 10"^ 


2.8 


5.5 


4.4 


3.3 


0.7 


1.0 


3.1 


0.93 


1.3 


2 D 


2.48 10^3 


6.3 


10.6 


9.1 


5.4 


2.1 


1.0 


4.9 


0.91 


1.5 


3 A 


3.82 10~^ 


2.0 


3.1 


2.8 


1.2 


1.1 


0.3 


0.4 


0.93 


1.1 


3 B 


5.86 1D~^ 


1.5 


2.4 


2.1 


1.3 


0.2 


0.4 


1.2 


0.92 


1.2 


3 B 


1.93 10^^ 


2.5 


4.9 


3.9 


2.9 


0.5 


0.6 


2.8 


0.93 


1.8 


3 D 


3.78 10"^ 


5.6 


9.1 


7.7 


4.9 


1.0 


1.8 


4.5 


0.91 


2.9 


4 A 


2.36 10~^ 


2.9 


4.2 


4.0 


1.4 


1.3 


0.1 


0.3 


0.92 


1.2 


4 B 


7.22 10"^ 


1.6 


2.5 


2.2 


1.2 


0.5 


0.4 


1.1 


0.91 


1.7 


4C 


3.91 10^^ 


2.1 


3.6 


3.0 


2.0 


0.5 


0.4 


1.9 


0.91 


0.9 


4 D 


7.01 10"-' 


4.8 


8.0 


6.9 


4.2 


1.2 


1.0 


3.9 


0.93 


3.4 


5 A 


2.91 10"^ 


8.7 


8.9 


8.8 


1.1 


0.6 


0.8 


0.5 


0.92 


3.2 


5 B 


4.50 10"^ 


2.3 


2.9 


2.7 


0.9 


0.8 


0.3 


0.5 


0.91 


3.2 


5C 


8.09 10"^ 


1.7 


2.5 


2.2 


1.1 


0.1 


0.2 


1.1 


0.91 


2.6 


5 D 


4.54 10^2 


2.1 


4.7 


4.2 


2.2 


0.7 


0.7 


1.9 


0.90 


2.9 


6 A 
6 B 
6C 
6 D 


1.80 10^1 


4.6 


7.6 


7.5 


1.5 


1.0 


0.3 


1.1 


0.91 


3.3 



Table 6: Normalised 2-jet cross sections as a function of and ^ together with their relative 
errors and hadronisation correction factors. Other details are given in the caption to tabled The 
bin nomenclature is defined in table |2l At high small ^ values are kinematically disfavoured 
or forbidden. 
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150 <Q^ < 200 GeV^ 


jet 


2-^ 


et 


3-jet 


1 a 


1 b 


1 c 


1 d 


1 a 


1 b 


1 c 


1 d 


1 


jet 


1 a 


100 


16 


5 


1 


59 


19 


2 





26 


1 b 


16 


100 


12 


2 


22 


72 


12 


1 


30 


1 c 


5 


12 


100 


8 





19 


77 


6 


19 


1 d 


1 


2 


8 


100 








16 


78 


6 


2-jet 


1 a 


59 


22 








100 











21 


1 b 


19 


72 


19 








100 








30 


1 c 


2 


12 


77 


16 








100 





16 


1 d 





1 


6 


78 











100 


4 


3-jet 


1 


26 


30 


19 


6 


21 


30 


16 


4 


100 




200 < Q2 < 270 GeV^ 


jet 


2-^ 


et 


3-jet 


2a 


2b 


2c 


2d 


2a 


2b 


2c 


2d 


2 


jet 


2a 


100 


16 


4 


1 


58 


19 


2 


1 


25 


2b 


16 


100 


13 


2 


22 


71 


13 


1 


29 


2c 


4 


13 


100 


9 





20 


76 


8 


20 


2d 


1 


2 


9 


100 








14 


75 


8 


2-jet 


2a 


58 


22 








100 











21 


2b 


19 


71 


20 








100 








28 


2c 


2 


13 


76 


14 








100 





17 


2d 


1 


1 


8 


75 











100 


6 


3-jet 


2 


25 


29 


20 


8 


21 


28 


17 


6 


100 




270 < Q2 < 400 GeV^ 


jet 


2-^ 


et 


3-jet 


3 a 


3b 


3c 


3 d 


3 a 


3b 


3 c 


3d 


3 


jet 


3a 


100 


16 


5 


1 


59 


19 


3 





27 


3b 


16 


100 


13 


2 


20 


71 


12 


1 


30 


3c 


5 


13 


100 


8 





20 


77 


7 


21 


3d 


1 


2 


8 


100 








16 


77 


8 


2-jet 


3a 


59 


20 








100 











21 


3b 


19 


71 


20 








100 








30 


3c 


3 


12 


77 


16 








100 





20 


3d 





1 


7 


77 











100 


6 


3-jet 


3 


27 


30 


21 


8 


21 


30 


20 


6 


100 



Table 7: The statistical correlation factors given in percent between different Pt, obs bins of 
different jet observables inside bins 1 to 3 as estimated from the data. The bin nomenclature 
is defined in table [21 



23 



400 < Q2 < 700 GeV^ 


jet 


2-^ 


et 


3-jet 


4a 


4b 


4c 


4d 


4a 


4b 


4c 


4d 


4 


J 61 


4a 


100 


15 


6 


1 


58 


20 


3 


1 


28 


4b 


15 


100 


13 


2 


19 


70 


14 


1 


31 


4c 


6 


13 


100 


9 





21 


76 


7 


23 


4d 


1 


2 


9 


100 








15 


78 


8 


2-jet 


4a 


58 


19 








100 











23 


4b 


20 


70 


21 








100 








30 


4c 


3 


14 


76 


15 








100 





20 


4d 


1 


1 


7 


78 











100 


7 


3-jet 


4 


28 


31 


23 


8 


23 


30 


20 


7 


100 



700 < Q2 < 5000 GeV^ 


jet 


2-^ 


et 


3-jet 


5a 


5b 


5c 


5d 


5 a 


5b 


5c 


5d 


5 


jet 


5 a 


100 


16 


5 


3 


57 


21 


5 


1 


28 


5b 


16 


100 


13 


2 


20 


68 


14 


1 


30 


5c 


5 


13 


100 


9 





20 


71 


9 


21 


5 d 


3 


2 


9 


100 








23 


72 


7 


2-jet 


5a 


57 


20 








100 











18 


5b 


21 


68 


20 








100 








29 


5c 


5 


14 


71 


23 








100 





19 


5d 


1 


1 


9 


72 











100 


7 


3-jet 


5 


28 


30 


21 


7 


18 


29 


19 


7 


100 




5000 < Q2 < 15000 GeV^ 


jet 


2-^ 


et 


3-jet 


6a 


6b 


6c 


6d 


6a 


6b 


6c 


6d 


6 


jet 


6a 


100 


16 


6 


3 


53 


22 


7 





28 


6b 


16 


100 


12 


3 


20 


62 


15 


2 


32 


6c 


6 


12 


100 


7 





21 


58 


9 


22 


6d 


3 


3 


7 


100 








22 


67 


14 


2-jet 


6a 


53 


20 








100 











19 


6b 


22 


62 


21 








100 








28 


6c 


7 


15 


58 


22 








100 





25 


6d 





2 


9 


67 











100 


13 


3-jet 


6 


28 


32 


22 


14 


19 


28 


25 


13 


100 



Table 8: The statistical correlation factors given in percent between different Pr.obs bins of 
different jet observables inside bins 4 to 6 as estimated from the data. The bin nomenclature 
is defined in table [2l 
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Measurement 


Cis[Mz) 


Uncertainty 




experimental 


theory 


PDF 




0.1195 


0.0010 


+0.0049 
-0.0036 


0.0018 


24.7/23 




0.1155 


0.0009 


+0.0042 
-0.0031 


0.0017 


30.4/23 


O-NC ' 


0.1172 


0.0013 


+0.0052 
-0.0031 


0.0009 


7.0/5 


"■jet ''■2-jet 0"3-jct 
0"NC ' CTNC ' 0"NC 


0.1168 


0.0007 


+0.0046 
-0.0030 


0.0016 


65.0/53 



Table 9: Values of as(M^) obtained from fits to the individual normalised inclusive jet, 2-jet 
and 3-jet cross sections and from a simultaneous fit to all of them. Fitted values are given with 
experimental, theoretical and PDF errors as well as with the normalised of the fit. 
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Figure 2: Distribution of the selected events (solid dots) shown as a function of selection vari- 
ables in an extended domain: the inelasticity y reconstructed with the electron-S method of 
2-jet events (a); the invariant mass of the two leading jets M12 (b); the transverse momentum 
ratio in the laboratory frame PT,h/ PT,e of 2-jet events (c); the rf^l^ of the inclusive jets (d). 
The data are compared with weighted MC simulations, DJANGOH (solid line) and RAPGAP 
(dashed line). Vertical dashed lines indicate the positions of kinematical cuts. 
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Figure 3: The nomialised inclusive jet (a), 2-jet (b) and 3-jet (c) cross sections in NC DIS 
measured as a function of Q^. The ratio of 3-jet to 2-jet cross sections is shown in (d). The 
measurements refer to the phase-space given in table [U The points are shown at the average 
value of Q"^ within each bin. For the inclusive jets the present data (solid dots) are compared 
to HERA-I published data [3], here shown corrected to the same phase space as the present 
measurement and shifted in for clarity (open dots). The inner error bars represent the statistic 
uncertainties. The outer error bars show the total experimental uncertainties, defined as the 
quadratic sum of the statistical and systematic uncertainties. The NLO QCD predictions, with 
parameters described in the section[3]and corrected for hadronisation effects are shown together 
with the theory uncertainties associated with the renormalisation and factorisation scales, the 
PDF and the hadronisation (grey band). The ratio R of data with respect to the NLO QCD 
prediction is shown in the lower part of each plot. 
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Figure 4: The normalised inclusive jet cross sections measured as a function of the jet transverse 
momentum in the Breit frame in regions of Q^. The points are shown at the average value 
of Pt within each bin. Other details are given in the caption to figure [3l 
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Figure 5: The normalised 2-jet cross sections measured as a function of the average transverse 
momentum of the two leading jets in the Breit frame (Py) in regions of Q^. The points are 
shown at the average value of (P^) within each bin. Other details are given in the caption to 
figured 
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Figure 6: The normalised 2-jet cross sections measured as a function of the proton momentum 
fraction ^ in regions of Q^. The points are shown at the average value of ^ within each bin. 
Other details are given in the caption to figure [3l 
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a, from Normalised 2-Jet Cross Section 
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Figure 7: The Osii^r = \/ (Q^ + Pt)I'^) values determined using the normalised inclusive jet 
cross sections (a) and the 2-jet cross sections (b), each measured in 24 bins of and Pt- 
The error bars denote the total experimental uncertainty of each data point. In each plot, the 
solid line shows the two loop solution of the renormalisation group equation, resulting from 
evolving the as{Mz) obtained from a simultaneous fit of all 24 measurements. The inner band 
denotes the experimental uncertainties and the outer band denotes the theoretical uncertainties 
associated with the renormalisation and factorisation scales, the PDF uncertainty and the model 
dependence of the hadronisation corrections. 
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Figure 8: The as{Mz) values determined using the normalised inclusive jet cross sections 
measured in 24 bins in Q'^ and Pt- Other details are given in the caption to figure |7l 
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Figure 9: The as(M^) values determined using the normalised 2-jet cross sections measured in 
24 bins in and (-Pr)- The solid line shows the two loop solution of the renormalisation group 
equation, as{Mz), obtained from a simultaneous fit of all 24 measurements of the normalised 
2-jet cross sections. Other details are given in the caption to figure|7l 
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Figure 10: The as{Q) values extracted by fitting the dependence of the normalised inclu- 
sive jet cross section in different regions of (a); as{Q) values extracted by fitting the (Px) 
dependence of the normalised 2-jet cross section in different regions of (b); as{Q) values 
extracted from the normalised 3-jet cross section in different regions of (c). In each case, 
the solid lines shows the two loop solution of the renormalisation group equation obtained by 
evolving the corresponding central value of the as{Mz). Other details are given in the caption 
to figure |7l 
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Figure 1 1 : The minimal of the fit of the NLO prediction with jj,r = Xr ■ a/ {Q"^ + Pt)/"^ to 
the normalised inclusive jet cross section as function of Xj- for 23 degrees of freedom. Vertical 
dashed lines represent the range where the renormalisation scale is varied in order to estimate 
the impact of missing orders beyond NLO, while the full line indicates the nominal choice of 
the scale with Xr = 1. 
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Figure 12: The values of as{iir) obtained by a simultaneous fit of all normalised jet cross 
sections in each bin. The solid line shows the two loop solution of the renormalisation group 
equation obtained by evolving the extracted from a simultaneous fit of 54 measurements of 
the normalised inclusive jet cross section as a function of and Pt, the normalised 2-jet cross 
section as function of and (Py) and the normalised 3-jet cross section as function of Q^. 
Other details are given in the caption to figure |7l 
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